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Relationships are obtained connecting the two-nucleon overlap function of the eigenstates in the 
(^4 — 2) particle system with the asymptotic behavior of the two-body density matrix for the ground 
state of the A-particle system. This makes it possible to calculate the two-body overlap functions, 
spectroscopic factors and separation energies on the basis of a realistic two-body density matrix. 
The procedure can be used in describing the (e, e'NN) and (7, NN) reactions where the two-body 
overlap functions are a key ingredient in the analysis. 
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I. INTRODUCTION 



In this paper we examine the two-particle overlap functions in interacting many-body systems and derive general 
relationships connecting them with the ground state two-body density matrix. The procedure is based on the 
asymptotic properties of the overlap functions in coordinate space, when the distance between two of the particles 
and the center-of-mass of the remaining ones becomes very large. This work can be considered as an extension of 
the analysis presented in where the asymptotic behavior of the one-body density matrix and the single-particle 
overlap functions was examined. Such an extension to the two-body sector is of considerable interest in view of 
present-day experimental possibilities in electromagnetically induced two-nucleon knock-out reactions like (e, e'NN) 
and (-y,NN). 

00 \ The effects of ground-state NN correlations on two-nucleon knock-out reactions have been intensively studied and 
discussed. This includes the cases of the (e, e'pp) and (e, e'pn) reactions on le O and 12 C ]o]-|i"o|] and their contributions 
(for 12 C ) to the semi-exclusive (e, e'p) reaction as well as of electroinduced yl| and of photonuclear processes 
JLsl p"9[ . It has been shown in j|,|7| that while for (e,e'pn) reactions the central short-ranged correlations (SRC) play 
only a marginal role, the two-proton knock-out in (e, e'pp) processes exhibits a sizeable sensitivity to the ground- 
state correlation effects. It was pointed out in [|lO| that the most promising extraction of SRC effect shows up in 
the longitudinal structure function which may be studied in the super-parallel kinematics. It was found || that 
the knock-out of So pp-pair dominates the 16 O (e, e'pp) reaction and this gives good perspectives for extracting 
information on SRC in nuclei from such processes. It has been shown also pd| that the two-hadron knock-out is 
a substantial contribution to the (e, e') reaction mechanism above the quasielastic peak. The correlation effects on 
4 .ffe (e, e'd) d reaction cross-section have been studied injnj]. The important role of the tensor correlations on the 
photonuclear cross-sections has been pointed out in [fl3|-|l5|. It was concluded in jljj] that combined analysis of (j,p) 
and (e, e'p) reactions, together with new data from (7, NN) and (e, e'NN) processes will lead to better understanding 
in electromagnetically induced knock-out and that namely the pp-channel is more sensitive to the SRC Jll|[l8| . 

The two-nucleon overlap functions and their properties are reviewed widely e.g. in |20| . They have been used in 
JT(| as one of the components of the charge-current density, containing the information on the nuclear structure in 
studies of two-nucleon knock-out and transitions to the low-lying discrete final states of the residual nucleus. The 
overlap functions used in || are constructed using phenomenological single-particle wave functions and a correlation 
function from pT|| . As mentioned in [pi a more sophisticated treatment should be given in principle on the basis of 
full calculations of the two-nucleon spectral function where both long- and short-range correlations are consistently 
considered and hence spectroscopic factors are automatically included. In |22| nrs t calculations of the two-nucleon 
spectral function of le O have been performed. Long-range correlations are treated by a Dressed RPA and SRC are 
included in the pair removal amplitude by adding defect functions obtained from solutions of the Bethe-Goldstone 
equation for the finite nucleus. 

The aim of the present paper is to study to what extent observables important for two-particle emission, such as 
the two-nucleon overlap function, can be derived on the basis of the nuclear ground state two-body density matrix. 
Our attempt is inspired by the possibility to extract one-nucleon overlap functions from the one-body density matrix 
of the nuclear ground state that has been proved in Jlj] and successfully applied for calculating the overlap functions, 
spectroscopic factors and separation energies associated with the bound states of the (A — 1) — particle system ||,f|,|]]. 
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II. DENSITY MATRICES AND OVERLAP FUNCTIONS 



The one-and the two-body density matrices are defined in coordinate space as: 

p {1) {x,x') = (^\a+{x)a(x')\¥ A *>) , (1) 

and 

pW( Xl ,x 2 ;x[,x' 2 ) = {¥ A ^\a+(x 1 )a + (x 2 ) a(x[) a(x' 2 )\^) (2) 

respectively, where "f^) is the antisymmetric A-fermion ground state normalized to unity and a + (x), a(x) are 
creation and annihilation operators at position x. The coordinate x includes both the spatial coordinate r and the 
appropriate spin and isospin coordinates. The matrices p*- 1 -* and p^ are trace-normalized to the number of particles 
and of pairs of particles: 

TrpW = / pW(x)dx = A, (3) 



Tr p m * / p^( Xl ,x 2 )d Xl dx 2 = ^4r^ ■ ( 4 ) 
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Since p^ and p^ are hermitian matrices their eigenstates ipa form a complete orthonormal set, in terms of which 
and p( 2 *> can be decomposed as 

oo 

^W^E^WWiV). ( 5 ) 



p^Hxux^x'^x'z) =£A£ 2 )vi 2) *(zi,z 2 )V4%i,4) • (6) 

a=l 

The eigenfunctions ipa and eigenvalues A^ are usually referred to as natural orbitals and natural occupation numbers 
p3[ . For fermionic systems the antisymmetry of the wave functions imposes the constraint < < 1. In analogy, 

(2) (2) 

the eigenfunctions ipa (xi,x 2 ) are called natural geminals and the associated real eigenvalues \ a - natural geminal 
occupation numbers ]2~i| ]. As a consequence of the antisymmetry of the nuclear ground state, the \^ obey the 
inequalities: 

< \f ] <(A-l)/2, for A odd 

< Af ] < A/2 . for A even (7) 

In uncorrelated (non-interacting) systems only the values Aq' = 0, 1 appear for both the one- and two-body 
natural occupation numbers. Whereas the presence of correlations always reduces the maximal one-body occupation 
number, this is not the case for the two-body occupation numbers. The upper bound in ([?]) is actually only reached 
for systems which are in a sense maximally correlated, as e.g. the occupation number of zero-coupled pairs in the 
seniority formalism in the limit of large shell degeneracy. 

Of more direct physical interest is the decomposition of the density matrices in terms of the overlap functions 
between the A particle ground state and the eigenstates of the (A — 1) and (A — 2) particle systems, since overlap 
functions can be probed in exclusive knock-out reactions. 

Inserting a complete set of (A — 1) eigenstates \a(A — 1)) into cq.(Q) one gets 

pWf^'l^^'W^M, (8) 

a 

where tfa\x) = {a(A — l)\a(x)\^(A)) is the one-body overlap function associated with the state \a(A — 1)). The 
spectroscopic factors are then defined by the norm 
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5« = (tp a \<p a ). (9) 

The two-nucleon overlap functions are defined as the overlap between the ground state of the target nucleus <f( A ) 
and a specific state VE , L C ' > of the residual nucleus (C = A — 2) p0| : 

SaOci.a*) = (*i C) |a(^i)«(^)|* (A) ) • (10) 
The two-particle spectroscopic factor is analogously defined by the norm 

si 2) - (<£ a \<S> a ) (11) 

and the two-body density matrix reads 

p ( - 2) (x 1 ,x 2 ;x' 1 ,x' 2 ) = Y,Ki x ux 2 )^ a (x' 1 ,x' 2 ) . (12) 

a 

As in the case of the single-particle spectroscopic factors where Sa^ < AmL jjj, one can find that Sa < AmL- 
Therefore, both one- and two-particle spectroscopic factors cannot exceed the corresponding maximal natural occu- 
pation numbers. 

III. RELATIONSHIPS BETWEEN OVERLAP FUNCTIONS AND DENSITY MATRICES 

It has been shown in jp that the one-body overlap functions associated with the bound states of the (A — 1) system 
can be expressed in terms of the ground state one-body density matrix of the A nucleon system. For example, in the 
case of a target nucleus with J™ = + , the lowest (no £7) bound state overlap function is determined by the asymptotic 
behavior (a — > 00) of the corresponding partial radial contribution ptj(r,r) of the one-body density matrix: 

tPn lj{r) = 79 -j—, r-j- , (13) 

C no ij exp(-fc no y a) /a 

where the constants C no ij and k no ij are completely determined by pij(r,r). In this way, both f no ij(f) and k no ij 
define the separation energy 

e r = E {A ~ 1] E {A) - ^ fc "°' 3 (14) 

and the spectroscopic factor S no ij = (<p no ij I fnolj)- The procedure gives also the next bound state overlap functions 
with the same multipolarity if they exist and its applicability has been demonstrated in Refs. 

Similar procedure but applied to the two-body density matrix is of significant physical interest for analyzing prop- 
erties of transfer reactions when two nucleons are removed from the target's ground state leaving the residual 
C — A — 2 system in a state 2 \ The procedure is possible again due to the particular asymptotic properties of 
the two-body overlap functions (|l0|). They satisfy the general equation |20| : 

~ V? - ^Vl + v 12 + 4 2 >) $ a (n,r a ) = a(r 1; r 2 ) , (15) 

where e [ I ] = - E ( A \ V12 is the internal di-nucleon interaction and o{ri,T2) is the nonlocal residual source 

term which contains the interaction between the two extra nucleons and the C-nucleus. 

The condition to obtain the asymptotic solution of eq.(|l^) for the two-body overlap function are considered in 



25||26|. In 1 25 1 the restrictions on the two-body interactions in the case of the system of three nonrelativistic 
spinless neutral particles have been formulated and the coordinate asymptotics of the discrete spectrum wave functions 
have been studied. In [ po|j26| the asymptotic behavior of wave functions and overlap functions have been investigated 
generalizing the amplitudes of Merkuriev |^5| and different types of overlap functions necessary in calculations of direct 
transfer amplitudes have been shown. In case when a cluster of two like nucleons (neutrons or protons) unbound to the 
rest of the system is transferred with a simultaneous transfer of both nucleons, as e.g., in reactions 18 O — 2n^ 16 O, 
or 16 — 2p — > 14 C, the following hyperspherical type of asymptotics is valid for the two-body overlap functions 
i0|j5|]26|: 
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$( r ,i Z )_iVear P |-y-^-t^iP + -r»J j ^ + , (16) 

where r = r x — r 2 and R = (rj + 1*2) /2 are the magnitudes of the relative and center of mass coordinates, respectively, 
m is the nucleon mass and E is the two-nucleon separation energy 

E = - E^ = -E%. (17) 

In the case of a target nucleus with J K = + the two-body overlap function ( |Io| ) can be written in the form: 

^AJMixii^) = ^{® a JLs{Yl,r 2 )®Xs{<Tl,02)}j M , (18) 
LS 

where 



XSAfg (0-1,0-2) = |xi(oi) ®Xi(cr 2 )| 



SA/ S 

= ( ^ TOs i ^w S2 |5M s ) Xi msi (oi)xi ms2 (o 2 ) (19) 



and $ a jz,s'Afi( r ii r 2) is the spatially dependent part of the overlap function. As suggested in M] it is possible to 
perform a decomposition into angular momenta 1 = l r and L# (L = 1 + Lr), corresponding to the relative and center 
of mass coordinates: 

<&aJSLM L (r, R) = >, $aJSLlL R (r,R) y Lfl (R)®r ( (?) \ T „ ■ (20) 

IL R 1 J LMi 

Then the two-body density matrix has the form: 

pW(x 1 ,x 2 ;x' 1 ,x' 2 )=J2T,Yl PjSLlLn (.r,R;r',R>) x 

JM LS IL R S'L'l'L' R 
L>S'VL> R 

Asul*r (^i > CT 2 ; ?, R) ^s^l'i'L' R (°i ) 02 ; ? , R') , (21) 

where the radial part of the density matrix is 

PjSLlLn (r,R;r',R')=J2KjSLiL R (r,R)<S> a js>LTL> R (r',R') (22) 

S'L'l'L' R a 

and the spin-angular function is: 

Ajl{ LR (a 1 ,a 2 ;r,R) = {{Y LR (R)®Y l (?)\ ® X s(<ti,ct 2 )} (23) 

In eqs. (|l8|)-(p2"|) a is the number of the state of the residual nucleus with a given total momentum J. 
Let us further consider the diagonal part of eq.(|22]): 



P { jLiL R (r,R;r',R') =J2KjSLiL R (r,R)$ a JSUL R (r',R'). (24) 

a 

Eq. ( |l6| ) implies that for large r' = a and i?' = b a single term ao (corresponding to the smallest two-nucleon separation 
energy) will dominate in the sum of the right hand side of eq. ( |24|) . As a consequence the two-body overlap function 
<& aa jSLiL R { r i R) can be expressed in terms of the two-body density matrix: 

$ ao jSLiL R (r,R) = — 2 j—rr- (25) 

PjLilJ^R^^) (26) 



A^ea;p{-fc^/(6 2 + ia 2 )} (&a+l a 2)- 



■5/2 
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Here the asymptotic form (|l6|) of the two-body overlap function is used with a constant k which defines the two-nucleon 
separation energy (|i~7|). 

Therefore, if the two-body density matrix p( 2 > is known, the constant k entering eq.(26) can be determined by the 

(2) 

slope of the partial radial contribution Pjq^il ( r i-R; a J^) at large a and b, while the unknown coefficients N can be 

(2) 

obtained from the asymptotic behavior of the spatially diagonal part PjsLiL R ( a i b', a, b). In this way, the constant k 
will give information about the two-nucleon separation (|l7|), and the norm of the overlap function ( p6| ) about the 
associated two-particle spectroscopic factor. The relationships obtained can be useful in analyzing the processes such 
as (e, e'NN) and (7, NN) which are studied widely nowadays. Of course, the method will be reliable when realistic 
density matrices are considered. 

In conclusion, the derivation of eq.(^6|) makes it possible to find a solution of the complicated problem about the 
relationship between the two-body overlap functions and the two-body density matrix. The use of the asymptotics 
of the overlap function in the particular case of two like nucleons unbound to the rest of the system (|l^) can help 
in analyzing the overlap functions for two-nucleon knock-out processes on the basis of correlated two-body density 
matrix for the target nucleus. In our future work, which is now in progress, the density matrices obtained in the 
Jastrow-type model |27|]2S| ] are used for actual calculations of the two-particle overlap functions. 

The authors thank Dr. C.Giusti for the helpful discussion and the Bulgarian National Science Foundation for the 
partial financial support [contract 4>- 527]. 
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